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We present eight types of spatial optical solitons which are possible in a model of a planar waveg- 
uide that includes a dual-channel trapping structure and competing (cubic-quintic) nonlinearity. 
Among the families of trapped beams are symmetric and antisymmetric solitons of "broad" and 
"narrow" types, composite states, built as combinations of broad and narrow beams with identi- 
cal or opposite signs ("unipolar" and "bipolar" states, respectively), and "single-sided" broad and 
narrow beams trapped, essentially, in a single channel. The stability of the families is investigated 
via eigenvalues of small perturbations, and is verified in direct simulations. Three species - narrow 
symmetric, broad antisymmetric, and unipolar composite states - are unstable to perturbations 
with real eigenvalues, while the other five families are stable. The unstable states do not decay, 
but, instead, spontaneously transform themselves into persistent breathers, which, in some cases, 
demonstrate dynamical symmetry breaking and chaotic internal oscillations. A noteworthy feature 
,.0 ■ is a stability exchange between the broad and narrow antisymmetric states: in the limit when the 

two channels merge into one, the former species becomes stable, while the latter one loses its stabil- 
I^H , ity. Different branches of the stationary states are linked by four bifurcations, which take different 

■ forms in the model with the strong and weak inter-channel coupling. 

PACS numbers: 42.65.Wi; 42.65.Tg; 42.79.Gn; 42.82.Et 
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i. introduction and the model 

Self-trapping of light beams in the form of spatial solitons in planar waveguides, which was demonstrated exper- 
imentally about two decades ago is one of fundamental effects in nonlinear optics . The variety of the spatial 
solitons may be greatly expanded if the waveguide is equipped with a multi-channel structure [|J , that can be created 
by a permanent transverse modulation of the refractive index, or induced, in a photorefractive crystal, by a pair of 
laser beams (with the ordinary polarization) illuminating the crystal in transverse directions, while the probe beam, 
• that creates the soliton(s), is launched in the extraordinary polarization 3]. Besides their interest to optical physics, 
manipulations of spatial solitons by means of multi-channel trapping structures have a vast potential for applications 
£SJ ■ to routing of data streams [2|. 

The transmission of the electromagnetic wave with local amplitude i\){z, x) along axis z in the multi-channel planar 
waveguide is modeled by the nonlinear Schrodinger (NLS) equation. In the normalized form, the equation is 
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iip z +ip xx -V(x)TP + 6n(\iP\ 2 )ip = 0, (1) 

where the second term accounts for the transverse diffraction, V(x) represents the transverse modulation of the 
refractive index that induces the multi-channel structure, and term ~ £n(|-0| 2 ) represents the optical nonlinearity, 
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Sn (\4>\ 2 ) = J ^2|V'| 2 : with ri2 > 0, corresponding to the ordinary Kerr effect. In the model including the Kerr term and 
simplest transverse modulation, V(x) — Vq cos (2ttx/L), solitons trapped in channels and a possibility of switching 
them between the channels were studied in Ref. [2]. It was found that the entire family of fundamental (single- 
humped) solitons, with any value of integral power, P = \ip(x)\ 2 dx, is stable. A family of stable double-humped 
bound states of fundamental solitons exists in this model too, provided that the power exceeds a certain threshold 
value which grows with the increase of V$. A possibility of the switching, realized as carrying the fundamental soliton 
over into an adjacent channel, was also shown in Ref. Q, under the action of a local kick ("hot spot"), represented 
by an additional term ~ 8(x — xq)S(z — Zo)ip in Eq. JTJ), where xq is a midpoint between the channels. 

The same model as the one introduced in Ref. [2j, but with z replaced by time t was later considered as an 
effectively one-dimensional Gross-Pitaevskii equation for the Bose-Einstein condensate (BEC) trapped in the optical- 
lattice potential [H (the general topic of condensates trapped in optical lattices was reviewed in Refs. @). In the 
model with the sinusoidal transverse modulation replaced by that of the Kronig-Penney (KP) type, i.e., a periodic 
array of rectangular potential troughs, soliton states were studied in Ref. Q, and an effective nonlinear band structure 
in the same model was analyzed in detail in Refs. [f\. The KP structure is quite relevant for applications to optics, 
where it corresponds to the simplest step-index profile of the transverse modulation. 

In addition to employing the multichannel settings, the variety of spatial-soliton states can be greatly expanded by 
using media with competing self-focusing and self-defocusing nonlinearities. A well-known example of that is provided 
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by the cubic-quintic (CQ) nonlincarity. Optical nonlinearities of the CQ type were observed in chalcogenide glasses 
|8(, and in some organic materials [9(. Actually, the CQ nonlinear response of these media is induced by an intrinsic 
resonance , which also gives rise to nonlinear (two-photon) absorption [10]. Nevertheless, according to analysis reported 
in Ref . [ll| , effects of the loss may be neglected in physically relevant settings, as experiments in optical crystals are 
conducted over sufficiently short propagation distances (a few centimeters) [l], Q . It is also relevant to mention that 
the CQ nonlinearity was predicted [HJ and recently observed [HI in composite optical media (colloids). 
Equation (TTJ) with the CQ nonlinearity can be cast into the following normalized form [lil. [la]: 

ty, + 1>xx = ^(x)^-2|?M 2 V'+IV'lV (2) 
In the uniform medium, which corresponds to V(x) — 0, Eq. (T5|) gives rise to the well-known family of stable solitons 



(3) 



ip so i{x,z) = e lkz , _ 

\ 1 + y/i - 4fc/3cosh [2\/k. 

which is parameterized by propagation constant fc which takes values in interval 

< k < fc max = 3/4. (4) 

In an adjacent interval, 3/4 < k < 1, solutions exist too, in the form of continuous- wave (CW) states with a constant 
amplitude, 



^(z) = e lkz \Jl±VT^k (5) 
(the CW solutions exist for all k < 1). The total power of soliton ([3]) is 

V3. ( VZ + 2Vk \ 

^^-^vf-^k)- (6) 

CW solution ([5]) with the larger amplitude, V'cw' are stable, while solutions V'cw are subject to the modulational 
instability. At k = 1 , both branches of the CW solutions merge and disappear through the ordinary tangent bifurcation 
[l7j . Note that the width of soliton © diverges as 
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in the limit of k — » fc max = 3/4, see Eq. (fj|, and, accordingly, the soliton asymptotically goes over into the stable 
CW state, V'cw! m this limit. In terms of the bifurcation theory, the disappearance of the bright soliton at this point 
is explained by its merger with a dark-soliton solution to Eq. ([2]) with V = 0, which also exists at k < fc max = 3/4. 
Exactly at k = 3/4, the dark soliton degenerates into a front solution, which vanishes in one limit, at x — > ±00, and 
asymptotically coincides with the CW state in the other, 

A natural possibility, which is promising both for the exploration of fundamental properties of guided solitons and 
for potential applications, is to study localized states generated as a result of the interplay between the CQ nonlinearity 
and multichannel waveguiding structures. In Ref. [141 ] . solitons were investigated in the CQ model including a single 
guiding channel of a rectangular shape. A distinctive feature of the channel-trapped solitons in that model is their 
testability: while in interval ([4]), which hosts solution family ([3]) in the free space, the channel supports a single soliton 
state for each k, the full existence region of the solitons includes an additional interval, 3/4 < k < fc max (the value 
of fcmax depends on the depth and width of the guiding channel), where two different solitons are found for a given 
value of fc, both being stable (note that a guiding channel does not give rise to soliton bistability in the model with 
the cubic nonlinearity [Hj]). 

The combination of the CQ nonlinearity with a periodic array of guiding channels of the KP type was studied in 
Ref. 15]. In addition to single-humped (fundamental) solitons, many families of stable multi-humped (higher-order) 
solutions, which also feature the bistability, were found in that model. Independently, a model combining the CQ 
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nonlinearity and a periodic sinusoidal modulation function was introduced in Ref . , where similar results for soliton 
families were obtained. Note that both these versions of the model including the CQ terms and periodic potential 
function V(x) generate not only soliton families in the semi-infinite spectral gap, but also solitons in finite bandgaps, 
provided that the potential is strong enough. The finite-gap solitons are stable too, but they do not feature bistability. 

The most relevant setting for applications to the all-optical switching is based on two parallel guiding channels, 
rather than a periodic array. The respective model is based on Eq. ([2]) with the effective potential in the form of 

[0, \x\ < L/2 and M > D + L/2, . . 

1 j ~ \ -V , L/2 <\x\ <D + L/2, W 

where Vq, D and L are, respectively, the depth and width of each channel, and the thickness of the barrier which 
separates them. This model is also interesting in terms of the soliton dynamics, as it opens a way to study spontaneous 
symmetry breaking (SSB) in spatial solitons supported by the two-channel configuration. The SSB in settings based 
on double potential wells has recently drawn a great deal of interest in the studies of BEC, where it was realized 
experimentally [20| , and analyzed theoretically, using, chiefly, finite-mod e ap proximations, that replace the underlying 
Gross-Pitaevskii equation by a system of ordinary differential equations [21| . A conclusion produced by the analysis is 
that double- well systems with the self-attractive cubic nonlinearity (self-focusing, in terms of optics) give rise to stable 
asymmetric states (which is known as "macroscopic quantum self-trapping" ) through destabilization of symmetric 
states, while the systems with the self-repulsive (defocusing) nonlinearity feature destabilization of antisymmetric 
solutions. The latter mechanism also generates stable asymmetric states. 

Actually, a similar finite-mode analysis of the SSB of CW states in a nonlinear optical coupler (which, for this 
purpose, is essentially tantamount to the double-channel waveguide) was published much earlier in Ref. [22(. That 
work reported results for the model with the self-focusing or defocusing Kerr nonlinearity, as well as for a saturable 
nonlinearity. SSB in a double- well optical setting has been implemented in a photorefractive crystal (which corresponds 
to the saturable nonlinearity) (23|. 

The objective of the present work is to find symmetric, antisymmetric and asymmetric self-trapped states (i.e., 
spatial solitons, in terms of nonlinear optics) in the two-channel model with the CQ nonlinearity, based on Eqs. 
<j2j) and ([9|), and explore bifurcations linking those states to each other. An obvious difference from the previous 
works, both those relying upon the finite-mode approximation [2l| and works which analyzed the full two-channel 
configuration in the combination with the cubic nonlinearity, such as Ref. (23 |. is the competition between the 
self-focusing and self-defocusing nonlinearities, which drastically alters various states, their stability, and bifurcations. 

Note that dynamical switching of nonlinear localized beams in the same model as considered here was investigated 



in a recent work [25[, by means of a variational approximation and direct simulations. The analysis had revealed four 
different transmission regimes, depending on the total power of the beam. However, stationary states were not looked 
for in that work. 

The rest of the paper is organized as follows. In Sec. HH we present basic types of stationary spatial solitons which 
are possible in the model, and analyze their stability in Sec. IIII1 Depending on the symmetry and width of the 
solitons, we identify eight types of the states, five stable and three unstable (in the same configuration with the cubic 
nonlinearity, only three soliton species exist). The stability of each species is established in a rigorous form, via the 
computation of eigenvalues for modes of small perturbations around the solitons, and verified in direct simulations. 
For all unstable solutions, the instability eigenvalues are real. However, in direct simulations the instability does not 
destroy the trapped states, but rather transforms them into robust breathers, that may feature dynamical SSB, and 
chaotic intrinsic oscillations, in some cases. In Sec. IIV1 we demonstrate generic families (branches) of different soliton 
species, and bifurcations that link them together. The bifurcation diagrams are essentially different in cases of the 
strong and weak coupling between the two channels, i.e., roughly speaking, for small and large values of separation L 
between them, see Eq. @. The solution branches and bifurcation diagrams are displayed in two different forms, viz., 
as the total power versus the propagation constant, cf. Eq. and also as an effective asymmetry parameter of the 
state (it is defined below as per Eq. (fT3"|) ) versus the total power. The paper is concluded by Sec. [Vj 



II. BASIC TYPES OF SPATIAL SOLITONS IN THE TWO-CHANNEL SYSTEM 



In the model with the CQ nonlinearity, each waveguiding channel can support localized states of two different types, 
which we will call narrow and broad ones. The "narrow" state is nothing else but soliton ([3]) trapped in the channel, 
while its "broad" counterpart may be realized as a trapped fragment of CW state ([5]). In fact, the bistability of 
soliton states found in the CQ model with the single channel [l4| is accounted by the coexistence of stable solutions 
of these two types. In the two-channel system, these states may be combined in different ways. First, symmetric and 
antisymmetric spatial solitons, each built of narrow (Fig. [1} or broad (Fig. \2§ beams trapped in each channel, is easily 
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(a) (b) 



FIG. 1: (Color online) Typical examples of symmetric (a) and antisymmetric (b) spatial solitons built of narrow 
beams trapped in the two channels. These examples and those displayed below for other generic types of spatial 
solitons in the two-channel system are displayed for the two-channel potential with Vb = 3, D = 8, L = 1, and 
propagation constant k = 3.25, see Eq. (flOl) . Borders of the trapping channels are shown (here and in other figures) 

by vertical dashed lined. 
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FIG. 2: (Color online) Typical examples of symmetric (a) and antisymmetric (b) spatial solitons formed by broad 
pulses trapped in the two channels. Parameters are the same as in Fig. Q] 



found as a numerical solution of the equation obtained from Eq. (|2|) by the substitution of ijj(z,x) = exp (ikz) u(x) 
with real function u(x), 

u" = ku + V(x)u-2u 3 + u 5 . (10) 

Numerical solution of Eq. (|10[) was performed by means of the relaxation method based on the Newton-Raphson 
algorithm (formally speaking, this equation with V(x) taken in the form given by Eq. @ can be solved analytically, as 
explicit solutions in terms of elliptic functions are available in each region where V{x) takes a constant value; however, 
the continuity conditions for u{x) and u'(x) at points |a;| = L/2 and |i| = D + L/2 give rise to very cumbersome 
transcendental equations). 

As mentioned above, the NLS equation with the CQ nonlinearity in the free space, i.e., Eq. @ with V(x) = 0, 
admits not only exact solutions in the form of the bright solitons and CW states, as given by Eqs. ([3]) and (O, but 
also dark-soliton solutions. Accordingly, the antisymmetric bound states of two broad beams, a generic example of 
which is displayed in Fig. HJb), may be realized as a "curtailed" dark soliton (the one whose flat background was 
chopped off) trapped in the two-channel setting. 

Further, using a broad beam placed in one channel, and its narrow counterpart in the other, it is easy to find 
solutions to Eq. (|10[) in the form of asymmetric composite states, as shown in Fig. [31 The composite states, as well 
as the bistability of spatial solitons trapped in a single channel , are specific features of the model with the CQ 
nonlinearity (one may expect the same features in models with more general competing nonlinearities) , which are not 
supported by nonlinearities of non- competing types, such as those represented by cubic or saturable terms in the NLS 
equation. The signs of the two components of the composite spatial solitons may be identical or opposite, see Figs. 
[2[a) and [3](b). We will refer to these two species as unipolar and bipolar composite solitons, respectively. 

Asymmetric states that, as said above, are well known in models combining double-well potentials and cubic or 



saturable nonlinearity [2 11 ] — 24j, have their straightforward counterparts in the present model, in the form of a narrow 
or broad beam trapped in one channel, while the other channel is left almost empty. Typical examples of these 
(strongly asymmetric) states, which are different from the above-mentioned moderately asymmetric composite spatial 
solitons, are displayed in Fig. Below, we refer to them as "single-sided" states. 
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FIG. 3: (Color online) Typical examples of unipolar (a) and bipolar (b) asymmetric spatial solitons built as bound 
states of broad and narrow pulses trapped in the two channels. Parameters are the same as in Fig. [T] 




FIG. 4: (Color online) Generic examples of the single-sided states, with one channel carrying a narrow or broad 
pulse (they are shown, respectively, by solid and dashed curves in the left and right parts of the figure), while the 
other channel is nearly empty. Parameters are the same as in Fig. [T] 



A noteworthy feature of the single-sided broad-beam solution is that, for k approaching the maximum value, 
kmax = 3/4, above which soliton §5§ cannot exist in the free space, according to Eqs. (0| and (0, and is replaced by 
the spatially infinite CW state, V'cw ( see ^q- ©)' ^ ne broad-beam solution goes over not into something close to 
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but, instead, into a front state which is similar to exact solution ([8]), as shown in Fig. O Recall that, in the 



free space, the front solution exists at a single value of the propagation constant, k = 3/4, while an external potential 
may support a family of front solutions (26| (in Ref . (26| , fronts were represented by solutions to the Gross-Pitaevskii 
equation which included the cubic nonlinearity and a sinusoidal potential). In the present model, it can be shown 
that front states may also be sustained by a single potential well, in the absence of the second channel. We note that 
such a solution (which is stable, see below) was not reported in the previous analysis of the single-channel CQ model 

Examples of the eight different species of spatial solitons possible in the present model (symmetric and antisymmetric 
narrow and broad states, unipolar and bipolar composite ones, and broad or narrow single-sided solutions) were shown 
in Figs. [HQ] for a set of parameters at which the coupling between the channels is relatively strong (this is seen in the 
lack of strong separation between two components of the solitons in Figs. HJa), HJa), and[3](a)). In the limit of L — > 0, 
when the buffer layer that separates the channels disappears, and we thus return to the single-channel setting, the 
symmetric states continuously go over into single- humped solitons (narrow or broad ones), which were studied in Ref. 
\l4. In the same limit, the antisymmetric states gradually turn into a dipole (soliton-antisoliton pair) trapped in the 
remaining single channel. These transitions are illustrated by Fig. [51 which shows the total power of each symmetric 
and antisymmetric soliton family, P, versus L - from L = 0, that corresponds to a single channel of width 2D, up 
to large values of L, that correspond to isolated channels of width D. Note that the power of the narrow symmetric 
state in the limit of large L, i.e., the total power of the set of two fundamental solitons, is naturally very close to 
2P(L = 0), since the power at L = is that of a single soliton. On the other hand, it is natural too that powers of 
the broad states trapped in the single channel of the double width (at L = 0), and in two widely separated channels 
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FIG. 5: (Color online) An example of a stable single-sided broad beam, which is very close to the transition into a 
front-shaped state. Parameters of the two-channel potential are Vq = 3, D = 3, L = 6, and the propagation constant 
is extremely close to value k max = 3/4 which corresponds to the transition between the dark-soliton and CW states, 
via front configuration ([5]), in the free space (see Eq. (TJ|), k — 0.75 + 10~ 16 . 
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FIG. 6: (Color online) The total power of the antisymmetric and symmetric narrow (a) and broad (b) solitons as a 
function of thickness L of the buffer layer between the channels. Fixed parameters are Vq = 3, D = 8, and k = 3.1. 
The dashed vertical lines mark small values of L at which the stability is switched, see text. 

(at larger L) are nearly equal, as, for the nearly-CW configurations, P is proportional to the total width of the area 
filled by the wave. 

Unlike the symmetric and antisymmetric states, the evolution of asymmetric ones, i.e., the trapped beams of the 
unipolar and bipolar composite and single-sided types, with the decrease of L is not continuous. It was found that, 
at some small critical value of L, each asymmetric state performs a jump, either into a symmetric state (unipolar 
composite and single-sided spatial solitons), or into an antisymmetric one (the bipolar composite spatial soliton). 

Examples of antisymmetric states found in the limit of L = in the single channel are displayed in Fig. [7] Such 
dipole states were not considered in the framework of the previous work dealing with spatial solitons in the single- 
channel CQ model [ijj]. However, they are qualitatively similar to the so-called "dark-in-the-bright" solitons that were 
investigated as solutions to the Gross-Pitaevskii equation with a parabolic confining potential [27| • The peculiarity of 
the model with the CQ nonlinearity is that it provides for the coexistence of two different trapped dipoles at a common 
value of k, similar to the bistability of the trapped fundamental solitons. A noteworthy finding is that the stability 
of narrow and broad antisymmetric solitons trapped in a single channel is opposite to that of their counterparts in 
the two-channel setting, where the families of narrow and broad antisymmetric states are, respectively, stable and 
unstable (sec below): in the single channel, the narrow dipole is unstable, while the broad one is stable. In fact, the 
stability switch occurs at very small finite values of L, which are marked by dashed vertical lines in Fig. [5] 

In addition to the states trapped in the channels, we have also found solutions for beams trapped in the buffer 
layer between the channels. Such settings are known as nonlinear antiwaveguides [H, [HI]- We do not display those 
solutions here because, as may be naturally expected, they are subject to a strong instability. 



FIG. 7: (Color online) An example of antisymmetric spatial solitons (dipoles), narrow and broad ones (shown by the 
dashed and solid curves, respectively), trapped in the single channel of width 2D (edges of the channel are marked 
by dashed vertical lines). Parameters are Vq — 3, 2D — 16, and k = 3.1. The narrow dipole is unstable, while its 

broad counterpart is stable. 

III. STABILITY OF THE SPATIAL SOLITONS 

Stability of the eight families of spatial solitons found in the two-channel setting was investigated in two complemen- 
tary ways, viz., by numerical computation of stability eigenvalues from linearized equations for small perturbations, 
and by means of direct simulations of perturbed solitons. Both methods have produced identical conclusions concern- 
ing the stability or instability of each family of the stationary states. 

For the analysis of small perturbations, the solution was looked for as 

ip(x, z) = e ikz \u{x) + e (U(x) exp (az) + V*{x) exp (<r*z)] , (11) 

where u(x) is a real stationary solution pertaining to propagation constant k, the asterisk stands for the complex 
conjugation, e is a real infinitesimal amplitude of the perturbation with complex eigenmodes U(x) and V(x), that are 
associated to a (generally, complex) eigenvalue a. As usual, the stability condition is Re(cr) = 0, which must hold for 
all the eigenvalues. The substitution of expression l(TTj) in Eq. ([2|) and linearization yield the eigenvalue problem that 
can be written a matrix form, 

/ L + ta 2u\x) [1 - u 2 (x)] \f U \ , 

\2u 2 (x)[i-u 2 (x)] L-ia yv y / 

where the Sturm-Liouville operator is L = d 2 jdx 2 — V(x) + u 2 (x) [4 — 3it 2 (x)] — k. Numerical analysis of the stability 
spectrum was performed by the computation of eigenvalues of a large-size matrix approximating the operator on the 
left-hand side of Eq. (fT2]) . 

The following results have been obtained. Five species of the spatial solitons are stable in entire domains of 
their existence (i.e., all the respective stability eigenvalues have zero real parts, up to the accuracy of the numerical 
computation), with the exception of the above-mentioned stability switch of the antisymmetric states at very small 
values of L. Thus, the following families of stable solutions can be identified: 

(i) Narrow antisymmetric states, an example of which is displayed above in Fig. [TJb). 

(ii) Broad symmetric states, see an example in Fig. HJa). 

(iii) Bipolar composite states, see an example in Fig. [^b). 
(iv,v) Both narrow and broad single-sided states, see Fig. [4] 

Other three species, i.e., broad antisymmetric, narrow symmetric, and unipolar composite states, are unstable, also 
in their entire existence domains (except for the broad antisymmetric state at very small values of L, where it becomes 
stable, as mentioned above). Eigenvalues accounting for the instability of these three species are always real. In fact, 
the instability growth rate for the broad antisymmetric states (typically, it takes values a ~ 0.2) is usually larger than 
its counterpart for the narrow symmetric and unipolar composite states, with the same parameters of the model, by 
a factor ~ 10. 



8 



We note that the instability of the narrow symmetric states and the stability of their antisymmetric counterparts 
might be expected, as the general analysis developed in Refs. [2!| predicts that bound states of two solitons supported 
by an external potential are unstable if the phase shift between the solitons is Acj) = 0, and they may be stable if 
A(f> = 7r. This prediction does not apply to the symmetric and antisymmetric broad states, as they originate not from 
solitons proper, but rather from trapped segments of CW states, as discussed above. 

It is relevant to iterate that the stability of the antisymmetric states is swapped at very small values of L, which 
implies the transition from the two-channel setting to the single channel: as shown in Fig. [51 at L —> the broad 
antisymmetric solution becomes stable, while its narrow counterpart loses its stability. In fact, broad antisymmetric 
states trapped in the single wide channel resemble dark solitons, and their stability complies with the known fact that 
dark solitons are stable in the Gross-Pitaevskii equation combining the self-repulsive cubic nonlinearity and parabolic 
potential trap [30] . 

Predictions of the analysis based on the computation of the eigenvalues were checked in direct simulations of 
the evolution of perturbed spatial solitons, of all the eight species found above. The simulations of Eq. ([2]) with 
the respective initial conditions were performed by means of the usual split-step method. They identify stable and 
unstable species of the solitons in full agreement with the predictions of the linear-stability analysis. In particular, 
all the five species of the spatial solitons that are expected to be stable, listed above under numbers (i) through (v), 
indeed demonstrate very robust evolution in the presence of perturbations (examples of the stable evolution are not 
displayed here, as they do not reveal anything essentially new), except for the above-mentioned instability of the 
narrow antisymmetric state at L — > 0. 

Simulations of the perturbed evolution of the other three species, that should be unstable, reveal that the instability 
occurs indeed (with the exception of the above-mentioned stabilization of the broad antisymmetric state in the 
case of very small L, through the stability exchange with the narrow antisymmetric state). However, development 
of the instability does not destroy the solitons, but rather transforms them into persistent breathers, as shown in 
Figs. [81TTTI This outcome of the nonlinear evolution is observed despite the fact that the instability growth rate 
for small perturbations is always real, while the transition to breathers would be more straightforward in the case 
of an oscillatory instability, accounted for by pairs of complex conjugate eigenvalues. Accordingly, the period of 
the established oscillations of the breather depends, although weakly, on the size of the initial perturbation which 
destabilized the underlying stationary soliton (in the case shown in Fig. the period will be larger by a factor ~ 1.5 
if the amplitude of the initial perturbation is smaller by a factor of 10). 

The breathers formed from the unstable symmetric and antisymmetric solitons, such as those displayed in Figs. [5] 
and [51 feature not only the phase shift of it between the oscillations in the two channels, but also SSB (spontaneous 
symmetry breaking), as clearly manifested, in panels (b) of the figures, by the difference in the peak powers in the 
channels. Additional simulations demonstrate that, unlike the period of the oscillations, the size of the symmetry 
breaking, if measured as the difference between the maximum values of the peak powers, practically does not depend 
on the size of the initial perturbation. 

In most cases, the breathers formed from unstable solitons exhibit strictly periodic oscillations, as in Figs. [U(a) 
and 111! However, one species of the unstable stationary states, viz., broad antisymmetric ones, may transform itself 
into a breather which features (seemingly) chaotic oscillations. Some irregularity is seen in the evolution of the peak 
powers in Fig. [H^b). The increase of the coupling between the channels, i.e., a decrease in L, makes the oscillations 
manifestly chaotic, see Fig. [101 Note that the instability growth rate of the unstable solitons, <r, which give rise to 
chaotic breathers is much higher than in cases when periodic breathers are generated: for instance, it is, respectively, 
a 0.05 and a ~ 0.5 in Figs. [SlandfTUl On the other hand, Fig. fTUTb) suggests that the chaotic breather restores, 
on average, the originally broken dynamical symmetry between the two channels. Indeed, in that case the average 
values of the peak power in the channels are equal, both being 1.955. 



IV. BIFURCATION DIAGRAMS 



The global description of various stationary soliton families found in the present model, and of links between them is 
provided by diagrams that show the integral power of relevant solutions, P, as a function of the propagation constant, 
k, along with diagrams that display an effective asymmetry of the solutions versus P. The latter characteristic is 
defined as follows: 

e= p (13) 

[recall P = J ^ \ip(x)\ 2 dx + / + °° |?/>(2:)| 2 dx is the total power of the spatial soliton]. The form of the bifurcation 
diagrams turns out to be quite different in the cases of strong and weak coupling between the two channels, which 
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(a) (b) 

FIG. 8: (Color online) Transformation of an unstable narrow symmetric stationary soliton into a persistent 
breather, initiated by a weak initial perturbation: (a) the evolution of the local intentisty; (b) peak intensities in the 
two channels versus z. Parameters of the two-channel trapping potential are Vb = 3, D = 4, L = 2, and the 
propagation constant of the unperturbed state is k — 2.85. The calculation based on Eq. (fT2")l yields instability 

growth rate a « 0.05 for the underlying stationary soliton. 

correspond to small and large values of thickness L of the barrier between the channels. These two situations are 
presented below separately, with labels identifying the eight basic types of the solution branches in the plots as follows: 

A and B - symmetric stationary solutions of the broad and narrow types, respectively; 

A' and B' - antisymmetric solutions of the broad and narrow types, respectively; 

C and C - unipolar and bipolar composite states, respectively; 

D and E - single-sided broad and narrow solutions, respectively. 

A. Strong coupling 

Generic bifurcation diagrams for the system with a relatively strong interaction between the beams trapped in the 
two channels are displayed in Figs. 1121 and IT51 Following the usual convention, continuous and dashed curves in these 
plots represent stable and unstable solution branches. Discontinuities at turning points in panels [T27 a) and !13f a) are 
due to the (well-known) problem with poor convergence of numerical solutions in a vicinity of such points. For the 
same reason, the left, nearly vertical, segments of branch E in Figs. [T27 b) and fTSFb) (and also in Figs. [T4T b) and 
[TST b) below) were found with a low accuracy, therefore they are shown by dots. Note that values of k at the two right 
turning points in Figs. [UK a) and !13f a). as well as in Figs. [T4T a) and fTST a) below, are close but not equal. 

We stress that the bifurcations at points where branches A, B and C meet in Fig. [T27 a). as well as where A', B' and 
C' meet in Fig. I13f a). actually involve not three but four different branches (in compliance with general principles 
of the bifurcation theory [13 )• Indeed, Figs. [T2T b) andfTBTb) clearly demonstrate that branches C and C exist each 
in two copies, that are mirror images to each other: one is composed of a narrow soliton in the left channel and a 
broad one in the right channel, and vice versa. It is noteworthy that branches C and D, which represent, severally, 
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(a) (b) 

FIG. 9: (Color online) The same as in Fig. [5] (with the same values of the parameters) for an unstable broad 
antisymmetric soliton. In this case, the instability growth rate for the underlying stationary soliton is a sa 0.25. 

the unstable unipolar composite solitons and stable single-sided broad solitons, meet and disappear through a tangent 
(alias saddle-node [17[) bifurcation in Fig. fTOT a). Also worthy to note is the fact that branches D and E, which meet 
and disappear at the lower left turning point in Fig. fTCT a). are both stable, contrary to the "naive" expectation that 
one of them ought to be unstable (this fact was already discovered in the single-channel CQ model jl4|). 

The bifurcation diagram in the plane of P and e actually describes the SSB, in terms of stationary solutions. Note 
that the diagram in Fig. fT2Tb) features a loop, which is a characteristic feature of the SSB in models with saturable 
[22| and CQ nonlinearities. In the latter context, closed bifurcation loops were found in the discrete (lattice) version of 
the CQ model j3l|, and in a system of two linearly coupled NLS equations with the CQ nonlinearity [32j. The presence 
of the bifurcation loop implies that the broken symmetry is eventually restored, under the action of the quintic term. 
On the other hand, the diagram in the same plane for solution branches A', B' and C in Fig. [T37b) does not display 
a closed loop. In fact, two mutually symmetric branches C in this diagram terminate because of difficulties with 
the numerical continuation to larger values of P; on the other hand, the branch of broad antisymmetric solitons, 
A', can be continued easily, and the continuation does not reveal its stabilization at any particular value of P (in 
direct simulations, the broad antisymmetric soliton always turns into a persistent breather, cf. Fig. H^b)), hence it is 
plausible that the configuration observed in Fig. [TlTb) is not going to form a closed bifurcation loop. As concerns the 
continuation of branch A in Fig. [T27b) (and also in Fig. [Tilfb) below) to P — > 00, this happens with k — > fc max = 3/4, 

when state A becomes asymptotically similar to CW state V'cw' see ^q- ©• Similarly, state A' in Figs. [TSTbl. and 
also in !15f b) below, become similar, in the same limit to a dark soliton in the free space. 



B. Weak coupling 

At essentially larger values of the thickness of the buffer layer between the channels (L), which implies weak 
interaction between them, the character of the bifurcation diagrams changes qualitatively, as shown in Figs. 1141 and 
HT1 In this situation, the tangent bifurcation through which branches C and D disappeared in the case of the strong 
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(a) (b) 

FIG. 10: (Color online) The same as in Fig. [HI with a difference that the separation between the two channels is 
L = 2 instead of L = 1. The stronger coupling between the channels results in the formation of a chaotic breather. 
The instability growth rate for the underlying stationary soliton is cr « 0.5. 

inter-channel coupling does not take place, cf. Fig. [T27 a). Instead, there occurs another bifurcation that involves 
solution families C, C and D (in fact, family D of the stable broad single-sided solutions consists of two branches, 
one to the right of the bifurcation point, and the continuation to the left of it). This bifurcation manifests itself 
in the diagrams that involve the sets of both symmetric and antisymmetric states, i.e., (A, B, C) in Fig. IMf a), and 
(A',B',C) in Fig. [ijja). 

In the plane of (P, e), the continuation of branches D toward larger values of the power does not reveal any trend 
to forming a closed loop. Instead, the branches remain virtually identical to lines e = ±1, which is obvious in Figs, 
rjjjb) and fT5T b). This feature is explained by the propensity of the broad single-sided soliton to turn into the front 
state (which obviously has e = ±1), as shown above in Fig. [5] 

In both cases of the strong and weak coupling, the model gives rise to four bifurcations: those involving sets of 
branches (A,B,C), (A',B',C), (C,D) and (D,E) in Figs. QJandQS and sets (A,B,C), (A',B',C), (C,C',D) and 
(D, E) in Figs. [14]and[T5] The transition between the bifurcation diagrams corresponding to these two cases amounts 
to the well-known generic type of the rearrangement, when a branch splits off from a pitchfork bifurcation pattern, 
leaving behind a bifurcation of the tangent type (in the present situation, the detaching branch is C, as suggested 
by Fig. [TIT a)). We did not aim to find an exact value of L (for fixed Vo and D) at which the bifurcation changes its 
character. 



V. CONCLUSION 

The aim of this work was to explore spatial solitons that can be trapped in a planar waveguide which combines 
two fundamental ingredients, namely, the dual-channel configuration and competing (cubic-quintic) nonlinearity. We 
have found eight species of the solitons, including symmetric and antisymmetric ones of the broad and narrow types, 
unipolar and bipolar composite solitons, and the broad and narrow species of single-sided solitons. In contrast to 
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FIG. 11: (Color online) The same as in Fig. [8] for an unstable unipolar composite state. Parameters are Vq = 3, 
D = 4, L = 2, and k — 2.85. The instability growth rate for this state is a « 0.05. 



E D 




k Power 
(a) (b) 

FIG. 12: (Color online) Generic examples of the bifurcation diagrams in the strongly-coupled system, for Vo = 3, 
D = 8, L = 1. This set of the diagrams displays symmetric, unipolar composite, and single-sided branches of the 
spatial-soliton solutions, (a) The integral power (P) versus the propagation constant; (b) the asymmetry parameter 
(defined as per Eq. ([13)0 versus P. Labels A, B, C, D and E are defined in the text. 

this situation, previously studied models combining two-channel traps with non-competing nonlinearities (of the cubic 
and saturable types) gave rise to three species of the states (symmetric, antisymmetric, and single-sided asymmetric 
ones). The stability of the eight families was explored through the computation of perturbation eigenvalues and by 
means of direct simulations. It has been concluded that three families are unstable, viz., narrow symmetric, broad 
antisymmetric, and unipolar composite states, each of the other five species being stable as a whole. All instabilities 
against infinitesimal perturbations are accounted for by real eigenvalues. However, unstable solitons are not destroyed 
by the instability, but rather turn into robust breathers, which feature dynamical SSB (spontancousncss symmetry 
breaking), between the field intensities in the two channels. The breathers generated by strongly unstable broad 
antisymmetric solitons turn out to be chaotic (rather than periodic). It is interesting that the broad and narrow 
spatial solitons of the antisymmetric type switch their stability in the limit when the two-channel setting is going 
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(a) (b) 

FIG. 13: (Color online) The same as in Fig. [T31 but for branches of the antisymmetric and bipolar composite states. 

Labels A', B', and C are defined in the text. 




k Power 
(a) (b) 

FIG. 14: (Color online) The same as in Fig. [T2l but for the weakly-coupled system, with Vq — 3, D = 8, and L = 8. 

to merge into the single-channel one: the broad state becomes stable, while the narrow one is destabilized. Four 
bifurcations linking different branches of the stationary states have been found. The bifurcation diagrams may be of 
two different types, corresponding to the strong or weak coupling between the two channels. In the latter case, the 
bifurcation diagram involving the symmetric and unipolar composite states features a closed loop (as seen in Fig. 
[T2l) . An additional finding, which is relevant to the single-channel setting as well, is the stable front-shaped state (as 
shown in Fig. [5]). 

Numerous stable spatial solitons possible in this model may find applications to all-optical data-processing schemes. 
Moreover, the potential of the two-channel model is not exhausted by the analysis of the eight types of the trapped 
states reported in this work. Additional species of spatial soliton are possible in it, the form of symmetric, antisym- 
metric, composite and single-sided configurations based on stable broad dipoles (intrinsically antisymmetric states) 
trapped in each channel, as well as complexes involving a fundamental soliton in one channel and a broad dipole in 
the other. A challenging generalization is also possible for a two-dimensional model with the CQ nonlinearity and two 
parallel trapping channels. In that case, each channel may carry a fundamental two-dimensional soliton or a vortex. 




k Power 
(a) (b) 



FIG. 15: (Color online) The same as in Fig. [T3l with parameters identical to those in Fig. QjJ but also including 

branches D and E of the single-sided solitons. 
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